Kondo Lattice without Nozieres Exhaustion Effect 
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We discuss the properties of layered Anderson/Kondo lattices with metallic electrons confined in 
2D xy planes and local spins in insulating layers forming chains in z direction. Each spin in this 
model possesses its own 2D Kondo cloud, so that the Nozieres' exhaustion problem does not occur. 
The excitation spectrum of the model is gapless both in charge and spin sectors. The disordered 
phases and possible experimental realizations of the model are briefly discussed. 
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Famous exhaustion problem formulated by Nozieres Q 
is a stumbling stone on the way from exactly solvable 
Anderson or Kondo impurity model Q to the periodic 
3D Anderson lattice (AL) or Kondo lattice (KL) mod- 
els, which are believed to provide the generic Hamiltoni- 
ans for mixed valence and heavy fermion materials 0, . 
The problem arises already for concentrated Kondo al- 
loys where the number or localized spins iVj is compara- 
ble with the number of sites N = L n in the n-dimensional 
lattice. In this case the number of spin degrees of free- 
dom provided by conduction electrons in a KL is not 
enough for screening Ni localized spins. As an option 
a scenario of dynamical screening was proposed, 0, Q 
where only part of spins screened by Kondo clouds form 
magnetically inert singlets. The low-temperature state 
of such KL is a quantum liquid, where N s singlets are 
mixed with N — N s "bachelor" spins, which hop around 
and exchange with singlets thereby behaving as effective 
fermions. The Nozieres' exhaustion is measured by a pa- 
rameter pn — Ni/ pqTk (the number of spins per screen- 
ing electron). Here Tk ~ Pq 1 exp(— l/p^J) is the energy 
scale of Kondo effect, J is the exchange coupling constant 
in the single-impurity Kondo Hamiltonian, p is the den- 
sity of states on the Fermi level of metallic reservoir. 

Second obstacle, which does not allow the extrapola- 
tion of Kondo impurity scenario to KL is the indirect 
RKKY exchange Ijj> between the localized spins, which 
arises in the 2nd order in J or in the 4th order in V (hy- 
bridization parameter in the generic AL Hamiltonian). 
The corresponding energy scale is 
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where Xjf = N E q Xc(q) exp(iq • B. jf ) and Xc(q) is 
the spin susceptibility of the electron gas. The Fourier 
transform of Xc(q) is an oscillating function, which 
strongly depends on the distance Rjf. If I < at an 
average inter-impurity distance, and |/| ~ Tk, then the 
trend to inter-site antiferromagnetic (AF) coupling com- 
petes with the trend to the one-site Kondo singlet forma- 
tion (Doniach's dichotomy Q). 

This competition prompted several possibilities to by- 



pass the exhaustion limitations. According to a scenario 
offered in H, E|, in the critical region |/| ~ Tk of Do- 
niach's phase diagram, where the AF correlations are 
nearly suppressed by the on-site Kondo coupling, the spin 
liquid phase enters the game. This phase is characterized 
by the energy scale 



pol(T)=p J 2 X s H ,(T) 



where x 
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The condition 



, ,, is the spinon susceptibility. 
3 ' (Tk)>Tk is easily achieved both in 3D and 2D case. 
The Kondo screening is then quenched in the weak cou- 
pling regime at T>Tk, so that the spin degrees of free- 
dom remain decoupled from the electron Fermi-liquid ex- 
citations both at high temperatures T^>Tk and at low 
temperatures T<cTr- (Curie and Pauli limit for mag- 
netic response, respectively). At T^O the KL behaves as 
a two-component Fermi liquid with strongly interacting 
charged electrons and neutral spinons 10]. This scenario 
develops on the background of strong AF correlation. It 
includes the possibility of ordered magnetic phases with 
nearly screened magnetic moments and, in particular, the 
quantum phase transitions. Due to separation of spin and 
electron degrees of freedom, the Luttinger's theorem in 
its conventional Fermi-liquid form is invalid in this state: 
/-electrons represented by their spin degrees of freedom 
give no contribution in formation of the electron Fermi 
surface. Such state is referred as a "small Fermi surface 
regime" in current literature. 

Another scenario for small Fermi surface regime was 
proposed in [Tll | . This scenario appeals to systems where 
the magnetic order is either fragile or entirely absent due 
to magnetic frustrations (e.g., to triangular lattices). A 
spinon gap carrying unit flux of Z2 gauge field is expected 
to arise in spin subsystem, and this gap prevents forma- 
tion of Kondo singlets for a finite range of Tk ■ As a result 
the Nozieres' exhaustion does not occur, and fraction- 
alization of excitations into spin-fermions and electrons 
exists like in the previous case. A possibility of forming 
the spin liquid with (7(1) gauge group and spin density 
wave ground state has also been pointed out in [ill ]. 
In the present paper we propose another paradigm for 
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FIG. 1: (a) Layered lattice of spatially separated charges in 
planes and spins in chains, (b) A fragment of a chain with 
Kondo clouds formed as "shadows" in metallic layers. 

fermion fractionalization in Kondo lattices, which pos- 
sesses the generic properties of KL but is not subject to 
the exhaustion limitations. This paradigm may be re- 
alized in strongly anisotropic Kondo lattices, where the 
metallic electrons are confined in 2D planes interlaid by 
insulating layers containing magnetic ions. Then the 3D 
reservoir of screening electrons is defragmented into L 
planar reservoirs. Each plane still possesses the macro- 
scopic number of spin degrees of freedom ~ L 2 enough for 
Kondo screening, provided the concentration of magnetic 
centers per metallic plane remains small. The spin liquid 
features may be observed in these systems if the distri- 
bution of magnetic centers is also anisotropic, namely, if 
they form chains oriented in z direction, and the inter- 
chain interaction is negligibly small. 

Leaving the discussion of experimental realization of 
such systems for concluding section, we begin with the 
theoretical description of an ideal configuration, where 
all chains penetrate the stack in z direction (Fig^i). 
The AL Hamiltonian for the quasiperiodic model of con- 
duction electrons confined in metallic layers (xy plane), 
and magnetic ions localized in insulating layers between 
metallic planes is 

IWcr jo ^ ~ ' 

jl ko- 

Here k is a 2D wave vector, the discrete indices numerate 
metallic layers I with a lattice constant a|| and magnetic 
sites j along the chains with a spacing a z . The coupling 
constant V characterizes hybridization between itinerant 
2D electrons in a plane I and localized states in two adja- 
cent chains j, and j + 1 (Fig. lb). We treat the electrons 
in metallic planes in terms of Bloch waves Qko- , while the 
localized electrons are characterized by Wannier func- 
tions dj G . The lattice is quasiperiodic in a sense that 
the periodicity of magnetic sites in the xy plane is not 
demanded, but the average distance A between the impu- 
rities within a layer exceeds the radius of Kondo cloud, 
i.e. satisfies the condition A <C hvp/Tx (vf and Tr- are 
Fermi velocity of 2D electrons and Kondo temperature, 



respectively). There is no interaction between the chains 
under this condition, and a single chain represents the z 
component of excitation spectrum. On the other hand, 
all chains contribute to the ccy-component of the spin and 
charge response of the AL. The effects associated with the 
inter-chain exchange will be discussed in the concluding 
part. 

We came to a situation where L two-dimensional Fermi 
reservoirs, each with capacity L 2 , screen Ni magnetic mo- 
ments arranged in such a way that the effective concen- 
tration of these moments per metallic layer is rii = Ni jl? 
satisfies a condition n$ajj <C 1. This capacity is enough to 
form screening Kondo cloud for each magnetic site within 
a given layer I independently of all other sites belonging 
to the same layer. On the other hand, two magnetic 
ions localized one above another in neighboring insulat- 
ing layers j, j + 1 share the same metallic screen (see Fig. 
lb). Replicating these dimers along z axis, one we come 
to a system of spin chains, interacting with a system of 
metallic layers stacked up in the xy plane. Elimination of 
the hybridization term V in the Hamiltonian J3J in ac- 
cordance with the standard Schrieffer- Wolff procedure, 
results in effective exchange Hamiltonian for each chain, 

Ni 

H int = E/ ^kk'Si+l.kk' {§j + Sj+l) (4) 
j=l,k,k' 

(see Fig. Hb). Here s;+i, k k' = l4+i,ko-°< x < x ' Cl+1 ' k ' a ' > 
Sj = \d)j a cf aa idj^> , and an exchange integral is esti- 
mated as J - V£Vk>/U. 

Thus, the original model is reduced to the anisotropic 
KL formed by the system of ID spin chains penetrating 
the stack of 2D metallic layers. Each spin creates two 
Kondo clouds in adjacent planes, and two neighboring 
spins see each other through a metallic screen by means 
of indirect RKKY-like exchange. This exchange may be 
either ferromagnetic or antifcrromagnetic. In this work 
the latter case is considered in terms of the Doniach's 
dichotomy 0. In our model this dichotomy should be 
reformulated. Since the long-range AF ordering is im- 
possible in ID chain, two competing phases are Kondo 
singlet and spin liquid. Complete Kondo screening is not 
forbidden by Nozieres' exhaustion principle, since the 2D 
screening layer is available for each spin in the chain. The 
Kondo screening is characterized by the energy scale Tr . 
Thus, the competing phases in the anisotropic KL are the 
Kondo singlet phase and the homogeneous spin liquid of 
RVB type with the energy scale given by Eq. (J2J. 

In order to describe the Doniach-like phase diagram we 
adopt the method of 9]. Namely, we derive an effective 
action functional by integrating out all "fast" fermionic 
degrees of freedom with the energies ~ D , where 2D is 
the conduction bandwidth. The "slow" modes give us a 
hydrodynamic action. Due to strong quasi ID anisotropy 
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there is no need in appealing to the mean-field approx- 
imation. After elimination of conduction electrons with 
the energies Do > e > T in metallic layers, the cou- 
pling J is enhanced, J — ► J= l/pohiT/TV. The in- 
direct RKKY-like spin-spin interaction mediated by the 
in-plane electrons arises along the chains: 



TJ 
tl 



dd 



(5) 



Here I is defined in Eq. JIJ with Xjj+i — 
N^ 1 Xc(q) expi(q z a z ). Since V/f7 <C 1, we adopt the 
nearest-neighbor approximation for RKKY interaction. 

Up to this moment we treated the spin chains in a 
single-site approximation. This approximation is legiti- 
mate until T 3> Tk ~ /• To move further, we decouple 
the Euclidean action of the model (J3J , J^J 
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by means of the Hubbard-Stratonovich scheme |l2| in 
terms of the fields Aj,j±i ™ > J2a {d] a di±\ t(T + c.cj , 

<t>l Ek CT (4-i,k CT (^> + ^+i,-) + c-c) HEl. Here 

= d T — e(— iV) + /i and Pj~ c = d T — in/ (2/3) are bare 
inverse single particle Green's Functions (GF) for con- 
duction electrons and local spins, respectively, (3 = l/T. 
The field </> describes the single-site Kondo screening and 
the field A stands for the spinon propagation along the 
ID spin chain with AFM coupling. The single occu- 
pancy constraint dLd^f+dt, dj i=l is preserved at each 
site in the chain by the semi-fermionic transformation 
[lij . These two fields resolve the Doniach's dichotomy, 
because the long-range AF order is absent in ID. 

We appeal to the uniform resonance valence bond 
(RVB) spin liquid state [13| and treat the spinon modes 
as fluctuations around the homogeneous solution in a nn- 
approximation, Aj,j±i — 



A with 



£ 2 (/?)=/r 



A(r)A(-r)dr. 



For this sake we add and subtract A in the inverse 
GF. The non-local inverse spinon GF T)^ 1 = d T — 
Ajj±i — iirT/2. has to be expanded in terms of A — 
A. Now the two interacting components of bose-like 
modes in two-sublattice chain are spinons and Kondo 
clouds represented in effective action by A^+i j 

and <j)i,i+m<pi+m,i { m — 0,1), respectively. The charged 
0-mode acquires dispersion due to the non-locality of Jij , 
while the in-plane dispersion of conduction electrons in 
Kondo clouds is integrated out. The neutral spinon mode 
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FIG. 2: Loop expansion for non-local action @ Solid and 
dashed lines in IT2 and II4 stand for electron and spinon prop- 
agators, respectively. 



is dispersive by its origin. The action in these terms is 



A 



*ff - 2^ 



b U i(uj n )\ 2 |A ij+ i(w n )| : 



(7) 



Tr log^ 1 ) + Tr log (V- 1 (A jJ+l ) + G Q 4>*iAu±i + cc) 

As usual in ID systems, the spin and charge sectors in 
A e ff are separated. 

The last term in Q may be represented as a loop ex- 
pansion. Two first diagrams are shown in Fig. 2. To 
calculate the diagrams, we use the non-local spinon GF 
V'j j +r (uj n ) = (T>i oc — A) -1 with cosine-like dispersion 



exp 



In 



iV^n + A 2 



Here r numerates sites in the chain, oj n = 2ttT(ti + 1/4) 
on the imaginary axis 0]. IS characterized by two 
branch cuts at [A, +00) and (—00,— A]. In the limit 
A -C ttT it rapidly falls down with growing |r| as 
X>° i+r (w„) AH/(_j Wn )M+i. Thus the main con- 
tribution comes from T>j j(u> n ) = 



i/y/u 



A 2 and 

V% ±1 K) = (uJn/ v^ITa 2 " - 1)/A. 

The polynomial effective action after the loop expan- 
sion is performed is given by 




A| 



+Tr\og(G 1 ) + TrlogKV )- 1 } + 0(|</>| 4 ). 



(8) 



The polarization loops H2 and U4 are shown in Fig. 2. 
The action © is gauge invariant in accordance with 
Elitzur theorem , and spin and charge modes are sep- 
arated both in the 3D lattice and in the Fock space. 

To estimate A, we refer to the properties of spin chains 
with AF coupling ^J. The quasi- long-range-order in 
these chains may be treated in terms of boson excitations 
in Luttinger liquid (LL) or fermion pairs in spin liquid. 
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FIG. 3: (a) Disordered anisotropic Kondo lattices ; (b) For- 
mation of spin ladder from interacting chains. 

The spin susceptibility of a chain, (A + A~) w=0 ~ A 2 , 
acquires Pauli form at T* ~ 8J 2 /E F so we assume 
A ~ T* in our estimates. This means that even in the 
critical region of Doniach's diagram, Tk ~ T 2 jEp, the 
spins are "molten" into spin liquid at T ~ Tk, and the 
Kondo screening is irrelevant at low T . 

Evaluation of 112,114 in the limit 7rT^A gives 
nf ~p ln(A/T) and Ilf ~p /A 2 . This leads to reduction 
of indirect exchange, /=J[l+I/Am(A/TRr)] _1 . Since the 
Kondo screening is quenched at T*>Tx, we remain in 
a weak Kondo screening regime even in the LL limit 
at T— >0. The screening results in reduction of LL 
sound velocity, hv=Ia z . As to the in-plane charge ex- 
citations, the formation of Kondo clouds is quenched 
at T^>Tk, so instead of coherent Fermi liquid regime, 
{4> + 4>~)u)^o behaves as a relaxation mode ~ [—iui/T + 
aq 2 + ln.(A/Tjr)] , where F, a are numerical constants. 

These features of two-component electron/spin liq- 
uid manifest themselves in thermodynamics. The log- 
arithmic corrections ~ In - (T*/T) are expected in 
low-T Pauli-like susceptibility of isotropic spin chains, 
whereas the logarithmic corrections to the susceptibility 
of charged layers are quenched as InA/Tjf. The over- 
damped relaxation mode should be seen as a quasielastic 
peak in \o- The ID spinons contribute to the linear-T 
term in specific heat thus mimicking the heavy-fermion 
behavior, while the contribution of Kondo clouds is frozen 
at low T. 

In real anisotropic crystals one may expect formation 
of distorted and dangling chains (Fig. DJi) instead of an 
"ideal" lattice (Fig. 1). Distortion means shift of two 
neighboring Kondo " shadows" in a stack. This effect may 
be modelled by a random overlap factor Wj in RKKY 
integrals, Ij — Wjl. The dangling bond effect means 
Wj = 0. Bond disorder may be treated in terms of ran- 
dom AFM chains ^t|- According to this theory the dis- 
order results in transformation of singlet RVB liquid into 
a random-singlet RVB phase with arbitrarily long singlet 
bonds. In chains with broken bonds the gaps arise due to 
the finite length effect, so the short chain segments does 
not contribute into the low-T thermodynamics. With 
increasing impurity concentration, the Kondo clouds be- 



gin to overlap and two-leg ladders with diagonal bonds 
arise along with isolated chains (Fig. OJd). The Nozieres' 
exhaustion is still not actual for these clusters. With fur- 
ther increase of the concentration of magnetic sites, the 
Doniach's problem restores in its full glory. In case of FM 
coupling /, true long-range order emerges in spin chains, 
but the Nozieres' exhaustion is still quenched. 

One may point out the class of layered conduct- 
ing/magnetic hybrid molecular solids as an object for ap- 
plication of above theory. These crystals are formed by 
alternating metallic cationic layers and insulating mag- 
netic anionic layers with radicals [N(CN)2] _ as building 
blocks and transition metal ions as carriers of localized 
spins 18]. Organic cations with magnetic ions in such 
systems form ordered stacks. The problem is in preparing 
metallic layers with large enough Fermi surface to make 
Kondo screening effective and to find insulating networks 
with large enough distance between magnetic ions. It is 
worth noting that the above mentioned dicyanamide rad- 
icals with Mn ions form planar Kagome sublattice, thus 
being a promising object for realization of fractionalized 
Fermi liquid scenario proposed in . 
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